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Abstract 

This  paper  proposes  the  joint  space  time  signal  detection  algorithms  for  downlink  DS-CDMA  systems 
with  multiple-input-multiple-output  (MIMO)  multipath  wireless  channels.  The  proposed  algorithm  first 
decomposes  the  multipath  channel  into  a  set  of  independent  parallel  additive  white  Gaussian  noise  (AWGN) 
subchannels  each  with  the  signal  of  the  user  of  interest.  Signal  detection  can  then  be  achieved  on  a  path- 
wise-combining  basis.  According  to  different  multiple  access  interference  (MAI)  suppression  techniques, 
two  algorithms  respectively  referred  to  as  the  JST-MVDR  and  JST-COP  are  proposed. 
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I.  Introduction 

Multiple-input-multiple-output  (MIMO)  systems  using  antenna  arrays  have  received  significant  attention 
recently  because  of  their  capability  of  system  transmission  rate  improvement  in  rich  scattering  environment 
[1],  Diversity  techniques  are  often  adopted  to  enhance  the  channel  reliability  in  wireless  MIMO  systems.  To 
obtain  a  set  of  diversity  channels,  conventional  spatial  diversity  techniques  generally  assume  that  the  signal 
fadings  at  the  receive  antenna  elements  are  independent  and  the  optimum  diversity  can  be  achieved  by 
the  elementwise  maximum  ratio  combining.  However,  in  real  propagation  environments,  especially  when 
the  receive  antenna  elements  are  spaced  insufficiently  far  apart,  fading  effects  are  much  more  uncorrelated 
among  multipaths  than  among  the  receive  antenna  elements  [2].  With  this  understanding,  we  design  the 
diversity  combiners  based  on  the  signals  received  from  each  individual  path  instead  of  combining  the 
signals  received  from  the  antenna  elements.  By  exploiting  the  path  structure,  we  develop  a  space  time 
path-wise  combiner  for  the  downlink  direct  sequence  code  division  multiple  access  (DS-CDMA)  systems 
with  MIMO  wireless  channels  using  antenna  arrays. 

II.  System  Model 

Consider  a  downlink  MIMO  DS-CDMA  system  with  Mp  transmit  antennas.  Each  antenna  simultaneously 
transmits  K  data  streams  at  rate  ^  (symbols/sec.),  respectively,  for  K  users.  The  baseband  signal  of  the 
mth  transmit  antenna  is  thus  given  by 

N-l  K 

zm  (t)  =  ^2  ^2  X]  dm'k  (*)  Sm’k  (n)  9(t~  iiN  +  n )  Tc) ,  (1) 

i  n= 0  k—1 

where  drnk  (i)  denotes  the  ith  transmit  data  symbol  for  the  kth  user  of  the  mth  transmit  antenna,  and 
sm,k  (0) ,  •  •  •  ,  smy  {N  —  1)  is  the  iV-chip  spreading  sequence,  N  is  the  spreading  gain,  Tc  =  is  the 
DS-CDMA  chip  period,  and  g  ( t )  is  the  normalized  L  chip  long  pulse-shaping  waveform.  Due  to  multipath 
propagation,  signals  experience  fading  effects,  delay  spreads  and  angle  spreads  when  they  go  through  a 
wireless  channel.  Assuming  totally  P  paths  arrive  at  the  receiver,  after  coherent  frequency  down-conversion, 
the  complex  baseband  equivalent  representation  of  the  received  signals  can  be  expressed  as 

Mt  P 

x  (*)  =  ^2  &ppp  (*) Zm  (*  -  tp) + n  w  *  (2) 

m= 1  p— 1 

where  ap  denotes  the  spatial  signature  of  the  pth  multipath  corresponding  to  the  response  vector  of  the 
receive  antenna  array  with  respect  to  the  direction  of  arrival  (DOA)  of  that  multipath.  f3p  (t) ,  a  zero  mean 
complex  Gaussian  process,  is  the  fading  amplitudes  of  the  pth  multipath,  {rp}  is  the  path  delays,  and 
n  (t)  is  AWGN.  Also  in  (2),  p  =  1,  •  •  •  ,  P,  and  m  =  1,  •  •  •  ,  Mt-  The  fading  amplitudes,  though  time 
varying,  are  assumed  constant  within  a  Q-symbol  period  of  channel  coherence  time.  Sampled  p  times  per 
chip  period  Tc,  the  Q  receive  symbols  can  be  expressed  as 
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Mt  K 

X(*)=[  x  (zT) ,  •  •  •  ,  x((z  +  Q-l)r- f)  ]  =  ^^ABGTSm,fc4,t«  +  N,  (3) 

m=l  fc=l 

where  A  =  [  aj ,  •  •  •  ,  a p  ]  denotes  the  Mr  x  P  receive  spatial  signature  matrix  with  Mr  denoting  the 

number  of  the  receive  antenna  elements,  B  =  diag  {/3}  with  /3  =  [  P1  •  •  •  f5P  ]  is  the  P  x  P  fading 
matrix,  G  =  [  g(ri),  •••  ,  g  (rp)]  is  the  Lp  x  P  temporal  signature  matrix  with  g  (tp)  denoting 
the  sampled  version  of  g(t  —  Tp).  The  spreading-code  matrix  for  the  mth  transmit  antenna  of  the  kth 
user,  S m}k,  is  a  Toeplitz  matrix  constituted  by  the  spreading  sequence  {sm,k  (n)}  and  appropriate  zeros  to 
represent  the  convolution  between  the  pulse-shaping  waveform  and  spreading  sequence.  The  noise  matrix 
N  is  again  assumed  to  be  AWGN.  Furthermore,  throughout  this  paper,  we  assume  the  data  symbol  <lm,k  (z) 
is  BPSK  differentially-encoded. 

The  spatial  signatures  and  temporal  signatures  can  be  accurately  estimated  by  many  channel  parameter 
estimation  algorithms  [3].  Before  exploring  the  proposed  algorithm,  we  assume  that  the  spatial  signatures 
and  temporal  signatures  of  the  received  signal  are  known  to  the  receiver.  In  addition,  without  loss  of 
generality  we  assume  that  the  first  user  of  each  transmit  antenna  is  the  user  of  interest  and  the  spreading 
sequences  {sm,i  (rz) are  also  known. 

III.  The  Proposed  approach 

Shown  in  Fig.  1  is  the  structure  of  the  proposed  joint  space  time  signal  detection  algorithm  which 
consists  of  a  space  time  multipath  channel  decoupler  and  a  path-wise  signal  combiner.  The  multipath 
channel  decoupler  is  used  to  decompose  a  multipath  channel  into  a  set  of  independent  subchannels  and 
then  the  path-wise  signal  combiner  takes  over  to  detect  the  signal  of  the  user  of  interest.  According  to 
different  MAI  suppression  criteria  used  by  the  decoupler,  two  algorithms  are  developed.  One  is  referred 
to  as  the  JST-complementary  orthogonal  projection  (JST-COP)  algorithm,  and  the  other  is  referred  to  as 
the  JST-minimum  variance  distortionless  response  (JST-MVDR)  algorithm. 

A.  The  Space  Time  Multipath  Channel  Decoupler 

The  space  time  channel  decoupler  decomposes  a  multipath  channel  according  to  the  space  time  signature 
of  each  multipath.  To  obtain  the  space  time  signatures  of  the  channel,  by  stacking  all  columns  of  X  (z)  in 
(3),  and  making  use  of  the  rules  uec{XYZ}  =  (ZT®X)  •  vec{ Y}  and  uec{X i  +  X2}  =  i>ec{Xi}  + 
uec{X2},  we  can  rewrite  the  receive  space-time  vector  as 


(z)  =  vec  (X  (z))  =  ^2  dm,k  ( i )  Q m,k(3  +  nJST, 


where 


=  [  SJC  fcgi®aK 


,  Sm,fcgP  ®aP  ]  -  [  qm.fc.l 


and  <\m,k,P  =  fcg p  <S>  ap  is  referred  to  as  the  space  time  signature  of  the  pth  path  contributed  from 

the  mth  transmit  antenna  for  the  kth  user.  Based  on  the  space  time  signatures  of  the  user  of  interest,  the 
space  time  channel  decoupler  of  the  JST-COP  algorithm  isolates  each  path  by  using  the  COP  matrix  of  an 
appropriate  self-exclusive  matrix.  To  retain  the  pth  path  of  the  channel  and  null  out  the  others,  we  define 
the  space  time  self-exclusive  matrix  as 


The  corresponding  COP  matrix  is  then  expressed  as 

ttCop  t  nCOP  fnCOpNP 

^  p,m  A  Kolp,m  y^p^m  )  ’ 

Furthermore,  in  addition  to  isolating  each  subchannel,  the  space  time  channel  decoupler  performs  to  help 
suppressing  the  MAIs  and  scalarizing  the  vector  signal  in  (4).  To  do  this,  in  conjunction  with  the  space 
time  signature  of  the  user  of  interest,  the  weight  vectors  of  the  JST-COP  space  time  channel  decoupler  are 
given  as 

wCOP  _  TTCOPr,  n\ 

™p,m  p,m  \'J 

and  thus  the  outputs  of  these  decouplers  are 

Up,m  (*)  =  (w“p)H  xJST  (z)  «  dm,  1  (*)  /3p  +  n£°p.  (8) 


By  collecting  all  up>m  as  um  ( i )  =  [  ui>m  (i) ,  «p,ra  (i)  ] T  =  dm,i  (i)  /3T+n^OP,  signal  detection 

can  be  performed  by  combining  the  signal  of  the  user  of  interest  through  the  differential  decoding  as 


dm,  1  (*)  =  -sign( (i  -  1)  u,„  (i)). 


Alternatively,  the  space  time  channel  decoupler  can  also  be  implemented  based  on  the  MVDR  filters.  The 
MVDR  filters  have  been  deployed  in  many  signal  processing  applications,  such  as  spatial  beamformers 
and  frequency  filters,  which  seek  to  find  an  optimum  weight  vector  so  as  to  minimize  both  the  interference 
and  the  noise  at  the  filter  outputs.  To  design  an  MVDR  space  time  channel  decoupler  for  our  purpose,  we 
define  the  corresponding  constraint  matrix  and  response  vector,  respectively,  as 


—  Qm, i,p  pf^p 


and 


f  =  1  0 


IPxl  ' 


The  MVDR  space  time  channel  decoupler  w?)  m  can  then  be  formulated  as 

arg  min  E  j  \upprn  (i)|2| ,  s.t.  Cpmw  =  f, 


v  p,m 


(9) 

(10) 

(11) 


where  up,m  (i)  =  w^mxJST  (i)  is  the  output  of  the  space  time  channel  decoupler  wpm.  We  may  readily 
solve  (11)  by  using  the  constraint  optimization  method  and  obtain 


=  (R,,)  1  CP 


iH 


(Ruj)-1  Cp 


(12) 


where  R„,  =  E  |xJST  ( i )  xJST  (i)H  |  is  the  space  time  covariance  matrix  of  xJST  (*).  The  output  of  the 
MVDR  filters  can  thus  be  expressed  as 


R p,m  (^)  ~  fdpdm,  1  (t)  +  (13) 

where  np^m  =  w;^m  nJST  (i)  is  the  noise  at  the  MVDR  filter  output.  From  (8)  and  (13),  it  is  obviously 
that  the  space  time  channel  decoupler  successively  decomposes  a  multipath  multiuser  channels  into  a 
set  of  independent  AWGN  subchannels.  Combiner  with  differential  decoding  can  again  extract  the  phase 
difference  between  neighboring  dnhi's  to  eliminate  the  phase  ambiguity  caused  by  /3  . 

IV.  Simulations  and  Discussions 

Consider  a  Mr  =  2  and  Mr  =  10  system  with  half  wavelength  spaced  uniform  linear  array  (ULA)  at  both 
the  transmitter  and  the  receiver.  Differentially  encoded  data  sequences,  each  spreaded  with  a  gold  code 
sequence  of  31-bit  long,  N  =  31,  are  transmitted  over  the  channel.  We  assume  a  Rayleigh-fading  multipath 
channel  with  totally  3  paths  seen  at  the  receiver  for  the  user  of  interest.  The  direction  of  arrivals  and  delays  of 
the  three  paths  are,  respectively,  9  =  [  — 21.5 °  —20.1°  34°  ]  and  r  =  [  0.54TC  2.42TC  0.52TC  ]  , 

where  Tc  =  0.8 gs  is  the  chip  period.  The  average  fadings  of  the  multipaths  are  assumed  equal  and 
normalized  to  0  dB.  Within  the  time  span  of  our  interest,  each  path  is  assumed  to  have  a  randomly 
selected  but  constant  fading.  Five  hundred  Monte  Carol  trials  are  carried  out  as  the  number  of  multiusers 
are  set  at  K  =  1  and  K  =  5,  respectively.  Two  hundred  symbols  are  processed  in  each  trial. 
Fig.2-3compare  the  bit  error  rate  (BER)  of  the  proposed  algorithms  with  the  space  time  RAKE  receivers[4]. 
When  in  the  single  user  K  =  1  case,  since  there  is  no  MAIs,  Fig. 2  shows  that  the  proposed  JST-COP 
and  the  JST-MVDR  algorithms  perform  almost  the  same  as  the  JST-RAKE  receiver  while  significantly 
outperform  the  T-RAKE  receiver.  In  the  presence  of  MAIs,  with  significantly  better  performance  in 
interference  cancellation,  Fig.  3  shows  that  the  JST-MVDR  algorithm  outperforms  the  others  when  K  =  5. 


V.  Conclusions 

This  paper  proposes  a  space  time  algorithm  for  signal  detection  for  wireless  DS-CDMA  down-link 
MIMO  systems.  Signals  received  at  the  antenna  array  are  decomposed  in  the  space-time  domain,  and  then 
the  path-wise  combiner  is  applied  to  detect  the  data  of  the  user  of  interest.  The  JST-MVDR  algorithm 
generally  has  better  BERs  than  the  JST-COP  algorithm,  since  the  former  is  capable  of  suppressing  more 
interferences  including  ISIs  and  MAIs. 
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Fig.  1:  The  structure  of  the  proposed  algorithm 
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Fig.  2:  The  BER  of  the  proposed  algorithm  when  K=1 
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Fig.  3:  The  BER  of  the  proposed  algorithm  when  K=5. 


